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Abstract. We describe an essential improvement of our recent algorithm for computing cohomol- 
ogy of Lie (super)algebra based on partition of the whole cochain complex into minimal subcom- 
plexes. We replace the arithmetic of rational numbers or integers by a much cheaper arithmetic of 
a modular field and use the inequality between the dimensions of cohomology H over any modular 
field IF P = Z/pZ and over Q: dim_H"(F p ) > dimf/(Q). With this inequality we can, by computing 
over arbitrary IF P , quickly find the (usually rare) subcomplexes for which dim//(lFp) > and then 
carry out the full computation over Q within these subcomplexes. 

We also present the results of application of the corresponding C program to the Lie superalgebra 
of special vector fields preserving an "odd-symplectic" structure on the (2|2)-dimensional super- 
manifold. For this algebra, we found some new basis elements of the cohomology in the trivial 
module. 



1 Introduction 

Recently we proposed a new algorithm for computation of cohomology of a wide class of Lie superalgebras. 
This algorithm reduces the computation for the whole cochain complex to a number of smaller tasks 
within smaller subcomplexes. One can demonstrate that if T is the computation time for the whole 
complex, then partition of the complex into N subcomplexes reduces the computation time roughly to 
the value T/N 2 . Thus, the approach appeared to be efficient enough to cope with several difficult tasks 
in computing cohomology for particular Lie (super)algebras [1-5]. More detailed experiments with the C 
implementation of the algorithm, including profiling, reveal that arithmetic operations over Q take the 
main part of computation time (usually more than 90% for large tasks) . The same is true if Q is replaced 
by 2 (though computation becomes somewhat faster). 

A standard way to reduce negative influence of this bottleneck is to compute several modular images 
of the problem with subsequent restoring the result over 7L or Q by the Chinese remaindering or an 
algorithm for recovering a rational number from its modular residue [6,7]. Though, as is clear, the sum 
of sizes of modules used for constructing images can not be less than the size of maximum integer 
in the final result, the modular approach allows the intermediate swelling of coefficients to be avoided. 
Moreover the use of modular images is much more advantageous in the case of (co)homology computation 
than in the traditional problems of linear algebra. As we demonstrate further, the overwhelming part of 
computation can be accomplished using only one modular image. Recall that the Gauss elimination, 
the basic constituent of algorithms for computing (co) homology, over fF p has only cubic computational 
complexity in contrast to the exponential one over Q or 2. 

Note that the approach presented here can be applied not only for the Lie superalgebras, but in more 
general case of computation of homology or cohomology, especially when there is a practical method of 
splitting (co)chain complex into smaller subcomplexes. 

To demonstrate the power of the new algorithm and program, we present the results of computation 
of cohomology in the trivial module for the algebra SLe(2). This is an example (for n = 2) of the Lie 
superalgebra of special (divergence free) vector fields on the (n|n)-dimensional supermanifold preserving 
the odd version of symplectic (periplectic, as A. Weil called it) structure [8, 9]. These superalgebras, being 
super counterparts of the Lie algebras of Hamiltonian vector fields, are vital in the Batalin-Vilkovisky 
formalism, see [10]. 



2 Combining Splitting Algorithm with Modular Search 



The fcth cohomology is defined as the quotient group 

H k = Z k/ B k = Kcr d k/ lm d k-i 

for the cochain complex 

Q , c o d° j d"-2 ck _ 1 C k C k+i _ _ _ ^ 

Here, the C k are abelian groups of cochains, graded by the integer k (called dimension or degree); the d k 
are differentials (d k o d k ~ 1 = 0); the Z k = Ker d k and B k = Im c? fc_1 are the subgroups of cocycles and 
coboundaries, respectively (see [11] for details). In order to apply without restrictions the linear algebra 
algorithms, we assume that the groups of cochains are additive groups of certain linear spaces or modules 
and we shall use the corresponding terms in the subsequent text. 

2.1 Splitting Algorithm 

To compute the fcth cohomology, it suffices to consider the following part of (1): 

C k ~ x C k C k+1 . (2) 



First of all we split (2) using the Z-grading in the cochain spaces induced by the gradings in the Lie 
(super)algebra (and the module over this algebra) involved in the construction of the cochain spaces: 

f C k-l dk ~] (jk fjk+1 \ _ ^\ ( C k-l d g \ (jk _ d £^ C k+1 J 

^ ' geG ^ ' 

Here, G C Z is a grading subset. 

It appears that, as a rule, any subcomplcx in a given degree g can be split, in turn, into smaller 
subcomplexes: 

Here S 1 is a finite or infinite set of subcomplexes. 

Equation (3) means that the spaces C l g split into the direct sum of subspaces 

and the matrices of the linear mappings dg can be represented in the block-diagonal form 

d g = (J) d l,s- 

ses 

The construction of these subcomplexes is the central part of the splitting algorithm. 
Thus, the whole task reduces to a collection of easier tasks of computing 

ff£ s Key 4.. Im<, : . (4) 

As a basis of the cochain space C k , we choose the set of super skew-symmetric monomials of the form 

c(e n ,...,e lfc ;a Q ) = e n A • • • A e 4fc <g> a a . (5) 

Here, e, and a Q are basis elements of the algebra and module, respectively, and e* is the dual element to 
ej, that is, e*(ej) = 5*-. The degrees of factors in (5) satisfy the relation 

gr(e 11 ) + • • • + gr(e lk ) + gr(a Q ) = g. 
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Notice that gr(e l ) = — gr(ej) and this is a serious obstacle to extraction of finite-dimensional subcom- 
plexes for infinite-dimensional Lie (super) algebras when computing cohomology in the adjoint module 
(important in the deformation theory). We also assume that i\ < ■ ■ ■ <ik- 
To construct a subcomplex 

Cts 1 ^ C ls ^ Cft 1 (6) 
from the sum in right hand side of (3), we begin with choosing somehow an arbitrary starting monomial 
m g, start °f the form (5). There are various choices of the starting monomial and the time and space 
efficiency of computation depends on the choice. Having no better idea, we use at present the following 
three strategies: choice of a lexicographically minimal, or lexicographically maximal, or random monomial. 
We call these strategies bottom, top and random, respectively. Among these, the top strategy seems to 
be most efficient (see experimental data in Tables 1 and 2) and it is used by default. Nevertheless, other 
strategies help sometimes to push through difficult tasks when the top strategy fails. 

Then we construct the three sets M* 1 " 1 , M k s and M^ 1 of basis monomials for C*" 1 , Cg S and C^+S 
respectively, by the procedure ConstructSubcomplex presented on page 11. 

The function TakeMonomialFromSet called within ConstructSubcomplex takes the current 
monomial from a set of monomials. 

The function InverselmageMonomials generates the set of (q — l)-monomials whose images with 
respect to the mapping a"^ 1 contain a given g-monomial. 

The function ImageMonomials generates the set of (q + l)-monomials whose inverse images with 
respect to the d q g s contain a given (/-monomial. 

In the finite-dimensional case, the loops in the procedure ConstructSubcomplex are finite and we 
obtain in the end a minimal subcomplex of the form (6). This is the unique minimal subcomplex involving 
the starting monomial start . 



2.2 Modular Search 

Let us consider in more detail the procedure of computation of cohomology within the subcomplex in 
accordance with formula (4). From now on we assume that C*" 1 , C* s and C^t 1 in (6) are finite- 
dimensional spaces over Q or F p or modules over Z. 

Since important in mathematics and physics fields 1R and C are, in principle, non-algorithmic objects, 
our main interest will be focused on the cohomology over the field Q (or its algebraic extentions). In 
accordance with a general theorem in the homological algebra, called the universal coefficient theorem 
[12], (co)homology with coefficients from an arbitrary abelian group G can be expressed in terms of 
(co)homology with coefficients in Z. Thus, we can carry out the computation over Z and then go to the 
coefficient group we are interested in. Let us consider now the connection between Hjf S (Q) and H k S (Z). 
In the finite-dimensional case, the group H k S (JL) is a finitely generated abelian group having the following 
canonical representation 

free part torsion 

u k gs {iL) ~ ( ze---ez )ez tl e---ez tr . (7) 

Here, (3 k , the number of copies of the integer group Z, is called the rank of the abelian group Hg S (Z) or 
the Betti number. The cyclic groups Z ti are called the torsion subgroups and their orders t t , having the 
property tj > 1, ti\t 2 , t 2 \t 3 , . . . and so on, are called the torsion coefficients. 

In the case of cohomology, the universal coefficient theorem is expressed by the following split short 
exact sequence 

- Hl s {-K) H k JG) - Tor^+^Z), G) - 0, (8) 

where the operation Tor is the periodic product of abelian groups. In our context 

Tor^+^Z),^) = {torsion iJ^+^Z)) <8> {torsion G). 

The term "split", in application to sequence (8), means the possibility to construct the isomorphism 

H^ s (G) ~ Hl s (Z) ®G® Tor^t 1 (Z) , G) . (9) 

Replacing G by Q in (9) and taking into account that Tor(A, Q) = for any abelian group A, we have 
#£ S (Q) ^ H k s {W.) ® Q. Since Z m <g> Q = for arbitrary m and Z <g> Q ~ Q, the dimension of #£ S (Q), 
interpreted as vector space over Q, coincides with the rank (Betti number) f3 k of the group H k S (Z). 
Our modular approach is based on the following 
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Theorem 1. 



Remarks: 



dimff g fc s (F p )>dim#£ s (Q). 



(10) 



1. Inequality (10) means that non-trivial cohomology classes computed over the field of rational numbers 
Q can exist only in the subcomplexes with non-trivial cohomology classes computed over the finite 
field IFp with arbitrary prime p. 

2. H. Khudaverdian turned author's attention to the fact that inequality (10) can be deduced immedi- 
ately from the universal coefficient theorem: considering the product Hg g (^E.) (g) JF p and taking into 
account representation (7) and isomorphism Z ® W p ~ JF p , we see that the dimension of Hg s (W p ), 
as a vector space over W p , can not be less than (i k (only additional dimensions may appear, if the 
torsions in Hg S (K) or i?^+ 1 (2) contain cyclic groups of the form 7L p m). 

Nevertheless, we give here a direct constructive proof in order to demonstrate in parallel the main 
ideas of (co) homology computation. 

Proof. To prove inequality (10), we have to compute (4) in such a way as to avoid cancellations of integers 
and apply the modular homomorphism <f) p : 1L — > F p at the end of computation. Thus, it is convenient 
to consider (4) over 2 instead of Q. 

We assume that p is odd and use a symmetric representation of 1F P , i. e., 



F p = 



1 



.,-1,0,1, 
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We will also apply <j> p component-wise to multicomponent objects over Z, like vectors and matrices. 
We begin with the following setup: 

- Cg" 1 , Cg s and C^ 1 are represented as finite-dimensional modules M~ = Z™, M = Z m , M+ = 
1} ' , respectively, i. e., dimC^J 1 = n, dimC^ = m, dimC^ 1 = I. 

— the differentials d^" 1 and dg S are represented (in the monomial bases of the form (5), in our case) as 
integer m x n and I x m matrices 



D' = m D'~- 



(d')T 



(df n 



and D = l D= 



d\---dl 
d[ ■ ■ ■ d l m 



respectively. We write *A to indicate that matrix A has i rows and j columns. 
— the matrices D and D' satisfy the relation DD' = 0. 

The computation of cohomology, i.e., construction of quotient module, can be reduced to the construction 
of so-called (co)homology decomposition [13] based on the computation of the Smith normal forms [14] 
of the matrices representing differentials. 

First of all let us determine the cocycle submodule, i.e., Ker D C M, by reducing the matrix D to 
the integer Smith normal form S = UDV. The matrix S has the form 



rS 

r 

-rQ 



r O 

m— r 
-rQ 



(11) 



where, »6 is the i x j zero matrix, r = rank^Z?, rS= diag (si, ... , s r ) , Si, . . . , s r are positive integers 
called the invariant factors of D. These invariant factors have the property Sj|sj + i for all i. Note that 
there is connection between the invariant factors and the torsion coefficients from formula (7), namely, the 
prime divisors of the invariant factors are also divisors of some torsion coefficients. The transformation 

/ m 

matrices U — IJJ and V — my are unimodular integer matrices, i.e., det U = ±1, dct V = ±1. With 
such determinants, these matrices are invertible and their inverses U^ 1 and V^ 1 are obviously integer 
matrices too. 

Now we should consider the coboundary submodule 



Im D' C Ker D C M. 
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Combining the relation 



SV^D' = UDD' = 



with the structure of the matrix S (see formula (11)), we can reduce the matrix D' determining cobound- 
aries to the matrix D' acting in the submodule of cocycles: 



V- l D' 



rO 

n 

m—r £)t 



Computing the Smith normal form S' = U'D'V' for the reduced coboundary matrix we get 



S' = m-rS'= 



r'S' r' O 

r n— r 

m—r—r'Q m—r—r r 



where r' = rank^D' and r'S'— diag (s[, ... , s' r ,) . 

m — r 

We can extend the transformation matrix V = rn-rjj' acting in the submodule of cocycles to the 
transformation matrix acting in the whole module M : 



U' = mU'= 



m — r 
r O 



-r m — rJJ' 



Here, r I is the r x r identity matrix. Using the transformation matrices U' and V we can transform the 
initial (monomial in our case) basis e = (ei, . . . , e m ) in the module M into the basis a = (ai, . . . , a m ) = 
eV {U'Y making the cohomology decomposition explicit 



coboundarie 



cohomology 



M = (ai, . . . , a r ) ® (a r+ i, . . . , a r+r >) ® (a r+r > +1 , . . . , a m ) . 

" v ' S v ' 

non — cocycles cocycles 

In this decomposition we have 

Kcr d g,s = (ar+1, ■ • • ,a m ) 

and 

Im 1 = ( ar +^ ■ ■ -> a r+r'} ■ 

The formula for the dimension of cohomology (Betti number) follows from decomposition (12) 



(12) 



dim s (Q) = (3 k = m -r-r'. (13) 

Now let us consider how (13) changes under <p p . The image of (13) takes the form 

dim H* s (W p )=m p -r p -r> p . (14) 

Since 4> p is a ring homomorphism, we have for arbitrary unimodular matrix A with integer entries 

det^p(A) = </> p (det A) = <f> p (±l) = ±1, 

that is 4> p maps the above transformation matrices into invertible matrices. Hence the number of elements 
in the decomposition basis a remains unchanged, m p = m. On the other hand, the invariant factors 
s'i, . . . , s' r i and s\, . . . , s r of the matrices S' and S divisible by p vanish, hence rank <fr p (S) — r p < 
r and rank 4> P (S') = r ' < r' and inequality (10) is proved by comparing (13) and (14). □ 
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2.3 Implementation 

An algorithm based on the above ideas was implemented in the C language. The program called LieCo- 
homology Modular has the following structure: 

1. Input Lie (super) algebra A, module X over A, cohomology degree (dimension) k and grade g. A and 
X should be defined over (some algebraic extension of) Z or Q. 

2. Construct the full set Mg of fc-cochain monomials in grade g. 

3. Choose a prime p for searching subcomplexes with non-trivial cocycles by computing over JF p . 

4. Choose an element m k g € Mg (the starting monomial). 

5. Construct a minimal subcomplex s : C^" 1 Cjj s C^t 1 sucn that rrig G Cg S . 

6. Compute n = dim H^ S (W P ). 

7. If n > 0, then compute over Z or Q (or their extensions). We can use here the Chinese remain- 
dering or the rational recovery algorithm as more efficient procedures than direct computation over 
Z or Q. 

8. Delete all basis monomials of Cg from Mg. 

9. If Mg is empty then stop computation, otherwise go to Step 4. 

In the current implementation we obtain the relations determining cocycles and coboundaries (in fact, 
the rows of matrices of differentials) within the procedure ConstructSubcomplex. These relations are 
generated one by one as by-product of the functions InverselmageMonomials and ImageMonomi- 
als. To prevent unnecessary memory consumption, every newly arising relation is reduced modulo the 
system of relations existing to the moment and, if the result is not zero, the new relation is added to the 
system. Thus, we automatically have the matrices of differentials in the normal form just after completion 
of the procedure ConstructSubcomplex. This process is obviously equivalent to the Gauss elimination 
method, the most standard method for the computation of the Smith normal form of a matrix. 

In recent years, a number of new fast algorithms for the determination of Smith normal form have 
been elaborated [15,16]. These algorithms appear to be well suited to the (co)homology computation. 
It is worth to study the possibility to incorporate these algorithms in our implementation. We could, 
for example, remove the generation of relations from the functions InverselmageMonomials and 
ImageMonomials making them as fast as possible. Then, after construction of subcomplex with the 
help of these modified functions, we should generate the matrices of differentials separately and apply the 
fast algorithms to these matrices. Of course, this modification should be done if the total computation 
time decreases without substantial increase in the memory consumption. The works [15, 16] contain a 
detailed analysis of the properties of the sets of primes most appropriate for application of modular 
algorithms to a given matrix. 

Here we give only a few comments concerning the choice of prime p in our algorithm. These comments 
are based mainly on experiments with the program. 

From the practical point of view, we should use only primes p for which all operations in IFp can 
be done within one machine word. Thus, for 32bit architecture we should choose p from the set of 8951 
primes (3, 5, ... , 92681). A good choice should not produce excessive cocycles. Of course, such cocycles 
will be removed at Step 7 anyway, but at the expense of additional work. In our context, an "unlucky" 
prime is the one which divides the invariant factors of matrices of differentials (or, in other words, the 
torsion coefficients of cohomology over Z) and, as is clear, the probability for a given prime to be unlucky 
diminishes as the prime grows. On the other hand, there is an increase of time (in the examples we have 
computed, up to factor 2 or 3) and space expenditures with increase of p within the set (3,5,..., 92681), 
so it makes sense not to use too large primes for searching subcomplexes with potentially non-trivial 
cohomology. In our practice, we use, as a rule, a compromise: a prime near the half of 32bit word, namely, 
p = 65537 = 2 24 + 1, i. e., the 4th Format number. 

However, quite satisfactory results can be obtained even with much smaller primes, as is illustrated 
in Table 3. The symbols in the boxes of this table mean that (non-zero) n = 6xm Wp ^ H(R(2))g, 

whereas H(H(2))g means the cohomology with coefficients in the trivial module for the Lie algebra H(2) 
of Hamiltonian vector fields on the 2-dimensional symplectic manifold. We performed computation over 
all modular fields JF 3 through F17. As is seen in Table 3, the results for Wn fully coincide with those for 
Q for all computed grades g G [—2, . . . , 8] (and for all cohomology degrees k). The table also illustrates 
Theorem 1, i.e., all boxes containing non-zero dimensions for the field (Q contain also non-zero (> same 
for Q) dimensions for all fields F p considered. 
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In Tables 1 and 2 we present (considering both algebra and superalgebra cases) the running times for 
computation over W 17 of cohomology if£(H(2)) (for k = 7, 4 < g < 8) and i^(SLe(2)) (for k = 6, < 
.9 < 4). The columns presented in these tables are: the dimensions of cochain spaces, i. e., the sizes of 
matrices of differentials; the running times in seconds for the top strategy of the choice of the starting 
monomial and comparison of the bottom and top strategies. 

The times in both tables were obtained on a 1133MHz Pentium III PC with 512Mb. Note that the 
maximum memory consumption is near 46Mb and near 14Mb for the tasks in Table 1 and in Table 2, 
respectively. 



Table 1. Timing for # 9 (H(2), Fit), k = 7 



9 


dim Cg 1 


dim Cg 


dimC* 9 fc+1 


Ttop 


Tfjottom /Ttop 


4 


1580 


1128 


479 


< i 


2.0 


5 


3382 


2730 


1388 


4 


2.8 


6 


6734 


6132 


3606 


27 


3.3 


7 


12766 


12818 


8546 


214 


3.7 


8 


23074 


25488 


18963 


1128 


4.5 



Table 2. Timing for fl£(SLe(2), Fi 7 ), k = 6 



9 


dimC^" 1 


dimC* 


dimC 9 fc+1 


Ttop 


Tbottom 1 Ttop 





1867 


6605 


22119 


i 


2.2 


1 


3528 


12162 


39796 


4 


1.1 


2 


6546 


22102 


70817 


21 


4.7 


3 


11878 


39652 


124768 


87 


5.7 


4 


21073 


70110 


217696 


413 


6.0 



3 Computing if*(SLe(2)) 

In this section we present the results of application of the program LieCohomology Modular to the 
Lie superalgebra of special vector fields preserving periplectic structure on (2|2)-dimensional superspace. 
Pcriplectic supermanifolds with a fixed volume element play an important role in the geometrical formu- 
lation of the Batalin-Vilkovisky formalism [17, 18], an efficient method for quantizing gauge theories. 

Recall that a periplectic or an odd symplectic manifold is an (n|n)-dimensional supermanifold equipped 
with an odd symplectic structure, that is an odd non-degenerate closed 2-form. In an analog of Darboux 
coordinates [19], it takes the shape 

n 

w = J2 dx * Ade *- ( 15 ) 

i=l 

Here, x 1 , . . . , x™ and 6\, . . . , 9 n are even and odd (Grassmann) variables, respectively. The vector fields 
preserving 2-form (15) form a Lie superalgebra denoted by Le(n). The elements of Le(n) can be expressed 
in terms of generating functions (also called hamiltonians) and these generating functions generate a 
nontrivial central extension of Le(n) called the Buttin algebra and denoted by B(n). The bracket for 
arbitrary two hamiltonians / and g is called the Buttin bracket or antibracket or odd Poisson bracket and 
takes the form 

i— 1 x / 

Here, p(f) is parity of function /. 
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The odd symplcctic structure is a super version of the ordinary symplectic structure u> = Y^i=i ^ ^ 
dpi, where q 1 and pi are both even. In symplectic case there is an invariant volume form p u = uo 11 
(Liouville theorem) and all the vector fields preserving u> are automatically divergence free. Contrariwise, 
in periplectic case the volume is not preserved (see [9], for more geometric consideration of the subject 
see [20]) and one can impose the divergence- free condition additionally. 

Thus, we come to the special Buttin algebra SB(n). Its generating functions / satisfy the divergence 
free condition 

Af = 0, (17) 

where A is the odd Laplacian 



A : 



E 

»=i 



dx l d9, ' 



(18) 



This Z\-operator is, actually, the Fourier transform with respect to the odd variables of the usual de Rham 
differential (see [9]). Together with its homology, A plays the key role in the formulation of so-called 
Batalin-Vilkovisky "master equation". We can slightly reduce the special Buttin algebra by removing 
constants from generating functions, i.e., taking the quotient algebra modulo the center Z. The resulting 
algebra SLe(n) = SB(n)/Z is called the special Leites algebra. 

Since the above algebras are infinite-dimensional, in order to compute cohomology, we introduce a 
grading by prescribing grades to the variables {x 1 } and {9i} with subsequent extension of the grading to 
the polynomial functions of these variables. Most natural grading can be provided by setting gr(x 4 ) = 1 
and gr(#j) = —1. Here the word "natural" means that, with this grading, the algebra B(n) (and Le(n) = 
B(n)/Z) contains the inner grading element 



(19) 



guaranteeing that all non-trivial cohomology classes lie in the zero grade cochain subspaces (see [11]). 
Unfortunately, there is no good inner grading element 1 in the divergence free algebra SB(n) (and SLe(n)). 
That is why the computation of cohomology for these algebras is much more difficult task than for the 
algebras without divergence free condition. 

Let us now turn to SLe(2). Its basis elements up to grade 1 are 

grades 
-2 

-1 




Here, x, y and 6, -0 are even and odd variables of the (2|2)-dimensional superspace, respectively; Ei and Oi 
are even and odd basis elements of the Lie superalgebra, respectively. Notice that odd generating function 
corresponds to the even element of superalgebra and vice versa. This is a property of Lie superalgebras 
with antibrackets called parity shift. 



1 In a private communication I. Shchepochkina suggested a general view on the inner grading elements for the 
algebras with antibracket. She noticed that ^" =1 giX l 0i can be used as the inner grading element in B(n) or 
Le(n) as well as (19) assuming that gi are arbitrary integers and gr(a; 1 ) = gi and gr(0j) = — gi 
grading element can belong to a divergence-free subalgebras SB(n) or SLe(n) only if 5Z™ =1 gi 



This inner 
(in virtue of 

equation (17)), i. e., the set {gi} must contain non-zero integers with opposite signs. This leads to impossibility 
to construct finite-dimensional cochain subspaces of a given degree. 
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Let us present also the initial part of the multiplication table in SLe(2) (given are only non-zero 
brackets) 



[E2, E 5 

[E2, Eq 

[E3, E4 
[E 3 , E 5 




[E5 , E4 

[Eq, E4 

[E 5 , E e 



= —2E4, 
= E 6 , 
= 2E 6 , 



[Oi,0 7 
[E 5 ,O r 
[E e ,0 7 

[o u o 8 

[E 4 ,O s 
[E 5 ,0 8 



= -E 2 , 

= -o 7 , 

= -Os, 

= E3, 

= -o 7 , 

= 8 , 



Studying this multiplication table we can obtain some information about the structure of SLc(2). For 
example, there are some subalgebras important in the construction of representations of SLe(2): 

— commutative negative grade subalgebra A <0 = (Oi, E 2 , E 3 ) ; 

— semisimple zero grade subalgebra A = (£4, E 5 ,E 6 ) ~ so(3) ~ sl(2) ~ sp(2); 

— non-positive grade subalgebra A< = A <0 ^ A , a semidirect sum of the semisimple algebra and the 
commutative ideal. 

The results of computation of cohomology Hg (SLe(2)) are presented in Table 4. The boxes of this table 
contain either three numbers (with possible indication of non-trivial cohomology class) or right arrow. 
The three numbers from top downwards are 

— dimCg , the dimension of the whole space of fc-cochains in grade g; 

— the number of minimal subcomplexes C^J 1 Cg C^.l 1 constituting the whole subcomplex 
in accordance with formula (3); 

— max se s dimCg ;S , maximum dimension of (k, g)-cochain subspaces among all the minimal subcom- 



The right arrow — > means that all subsequent boxes in the row contain the same information. The reason 
for this is that all the relations defining cocycles and coboundaries for fc-cochains coincide with those for 
(fc — l)-cochains multiplied by the 1-cochain c{0\) = c(8ip). 

In our computation we found four genuine cohomology classes, i. e., generators of the cohomology 
ring, a, /3, 7,5. They are parenthesized in the table. Their multiplicative consequences are underlined. 
Notice that the cohomology ring contains nilpotents and zerodivisors: there are arbitrary powers of the 
cocycle a = c(6ip), but [3a — 0, -fa — and Sa 2 = 0. 



4 Concluding Remark 

When computing cohomology we start with the construction of the full set of (fc, g)-monomials. At the 
moment we do not see how to avoid this in deterministic algorithms. To represent the set of monomials, 
we need to allocate n = Ix dim Cg elements of memory representing basis elements of algebra and module 
(I — fc for the trivial and I — k + 1 for any non-trivial module). In our implementation we represent basis 
elements by two-byte integers. For the last box in column 9 (fc = 9) of Table 4 we have dim Cg = 648308 
and the set of monomials occupies near 12MB. The dimensions grow very rapidly, so, in fact, we are 
working on the brink of abilities of 32bit architecture and, even theoretically, we can add only a few rows 
to Table 4. But, as to arithmetical difficulties, we have some progress. 



plexes. 
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Algorithm: ConstructSubcomplex 
Input: start! starting (k, p)-monomial 

Output: Mg" 1 , M k s , Mg^ 1 , monomial bases of cochain spaces in the current 

subcomplex s: C k ~ x ^» C k iS ^4 Cg+\ such that m 9Start € Mg S 
Local: M*" 1 C M g fc -\ M k t3 C M 9 fc s , Mg^ 1 C M k +\ working subsets 
of currently "new" (not yet processed) monomials; 
Wg,! 1 . Wg tS , Wg^ 1 , working sets of monomials; 
rn-gTs 1 , m g,s, mgts 1 ' working monomials 
Initial setting: 
1: M*" 1 := M*" 1 := 
2: M 9 fe s :=M 9 fc 3 := {m* iStart } 
3: M^ 1 := M^ 1 := 
Loop over fc-monomials: 
4: while M k ^ s / do 

5: m' >s := TakeMonomialFromSet(M k a ) 

Supplement the set Mg" 1 : 
6: W^g.s 1 := InverseImageMonomials(d k ~ a 1 ,mJ iS ) 
7: Mg 1 " 1 : Mg fc " 1 U W*! 1 
8: M*" 1 := M*" 1 U W*" 1 

Supplement the set Mg+ X : 
9: Wgi 1 := ImageMonomials(d k a ,m k s ) 
10: Mg^ 1 := Mg^t 1 U W^t 1 
11: Mg^ 1 := Mg^ 1 U W^t 1 

Exclude processed monomial m h g >s : 
12: Mg fc s := M 9 fe s \ {m* ,} 

Loop over (fc + l)-monomials: 
13: while M^ 1 do 

14: m^ 1 := TakeMonomialFromSet(M k + 1 ) 

Supplement the set M k s : 
15: Wg,s := InverseImageMonomials(dg s,™^ 1 ) 

16: Mg fc ' s := M 9 fe s U W^ s 

17: Mg fe s := M 9 fe s U M/g fe s 

Exclude processed monomial : 
18: M k + X :=M k + 1 \{m k g + 1 } 

19: od 

Loop over (k — l)-monomials: 
20: while M^ 1 / do 

21: rag" 1 == TakeMonomialFromSet(Mg~ 1 ) 

Supplement the set M k s : 
22: Wg fe s := ImageMonomials(d k ^ s 1 , m^ 1 ) 

23: Mg fc ' s := M 9 fe s U Wg% 

24: M 9 fc s := M 9 fe ' s U W k s 

Exclude processed monomial m k ~ 1 : 
25: Mg^ 1 := M*" 1 \ {m^ 1 } 

26: od 
27: od 

28: return M k ~\ M k . s , M^ 1 
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Table 3. dim// g fc (H(2), R) for R = Q and F p , p = 3, 5, 7, 11, 13, 17; (fc,g) G [1, . . . , oo) x [-2, . . . , 8] 



g\k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


-2 




I3I5I7 
In I13 
{IitIq} 






I3I5I7 
I11I13 
{IitIq} 














-1 




2 3 


2 3 


2 5 


2325 


2 3 



















Is 


Is 


3 3 


3 3 


I3I5I7 
I11I13 
{IitIq} 










1 




2 5 


2 5 


2 3 2 7 


83 2 7 


63 












2 


Is 


4 3 


3 3 


1 3 3 5 


8 3 4 5 


H3I5 


4 3 


Is 


la 






3 




2 7 


2 3 2 7 


63 2 5 


10 3 2 5 


10 3 


8 3 2 5 


6 3 2 5 


2 3 






4 






5 3 1 5 


20 3 4 5 


31 3 3 5 


17 3 1 5 3 7 


2 3 4 5 4 7 


4 3 3 5 l7 


3 3 






5 


2 3 


63 


4 3 


83 2 7 
2n 


34 a 2 7 
2n 


42 3 2 5 


20 3 6 5 2 7 
2n 


10 3 4 5 2 7 
2n 


63 






6 


Is 


4 5 


3 3 3 5 


83 3 7 


19 3 4 7 


35 3 1 7 


45 3 3 5 


38 3 4 5 


17 3 1 5 


3 3 




7 




2 3 2n 


16 3 2 6 
2n 


48 3 10 5 

2l3 


64 a 14 5 

2l3 


36 3 6 5 2 7 


30 3 2 5 2 7 

2l3 


60 3 2 s 

2l3 


44 3 2 7 
2n 


8 3 2 7 
2n 




8 


3 3 


83 


83I5I7 


27 3 4 5 4 7 


87 3 3 5 3 7 


110 3 5 5 l7 


73 3 18 5 5 7 

In I13 
{IitIq} 


65 3 15 5 3 7 


56 3 3 5 


I83I5I7 
In I13 
{IitIq} 


1 3 
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Table 4. H*(SLe(2)) for (fc, g) € [1, . . . , oo) x [-2fc, . . . , -2fc + 20] 



3 + 2fc\fc 


1 


2 


3 


4 


5 


6 


7 


8 


9 


k > 9 





1 

1(a) 


1 

1 a 2 
1 


1 

la 3 
1 


1 

1 a 4 
1 


1 

la 5 
1 


1 

1 a 6 
1 


1 

la 7 
1 


1 

1 a 8 
1 


1 

1 a 9 
1 


1 

l« k 
1 


1 


2 
2 
1 




















2 


3 
3 
1 


4 
4 
1 




_ 




> 








> 


3 


6 
6 
1 


12 

8 
3 












_ 






4 


8 
8 
1 


23 

13 (/3) 
4 


26 
13 
4 
















5 


10 
10 
1 


44 
16 
6 


56 
18 
7 




— ► 












6 


12 
12 
1 


73 
22 

8 


118 
26 
15 


121 
26 
15 














7 


14 
14 
1 


116 
26 
10 


226 
34 
23 


246 
34 
23 














8 


16 
16 
1 


171 
30 
14 


414 
39 
36 


491 
42 
41 


492 
42 
41 












9 


18 
18 
1 


244 
34 
18 


718 
48 
52 


952 
52 
71 


970 
52 
71 












10 


20 
20 
1 


333 
38 
23 


1182 
54 
80 


1780 
65 
124 


1867 

65 (7) 
124 


1867 
65 
124 










11 


22 
22 
1 


444 
42 
28 


1870 
60 
119 


3204 
72 
197 


3528 
76 
197 


3534 
76 
197 








> 


12 


24 
24 
1 


575 
46 
34 


2858 
66 
176 


5584 
84 
311 


6546 
88 
358 


6605 
88 
358 




_ 






13 


26 
26 
1 


732 
50 
40 


4224 
72 
241 


9398 
92 
489 


11878 
104 
606 


12162 
104 
606 










14 


28 
28 
1 


913 

54 
47 


6082 
78 
330 


15343 
100 
787 


21073 
113 
1009 


22102 
118 
1009 


22119 
118 
1009 






> 


15 


30 
30 
1 


1124 

58 
54 


8552 
84 
434 


24348 
108 
1187 


36540 
130 
1578 


39652 
134 
1598 


39796 
134 
1598 


— 






16 


32 
32 
1 


1363 
62 
62 


11766 
90 
570 


37649 
116 
1776 


61884 
140 
2556 


70110 

153 (5) 
2802 


70817 
153 8a 
2802 


70817 
153 
2802 






17 


34 
34 
1 


1636 
66 
70 


15892 
96 
721 


56848 
124 
2528 


102466 
150 
4127 


122154 
164 
4610 


124766 
170 
4610 


124798 
170 
4610 






18 


36 
36 
1 


1941 
70 
79 


21114 
102 
912 


84034 
132 
3583 


165999 
160 

6594 


209566 
184 
7510 


217692 
188 
7510 


217972 
188 
7510 






19 


38 
38 
1 


2284 
74 
88 


27622 
108 
1127 


121790 
140 
4900 


263372 

170 
10005 


353762 

196 
11624 


376010 
210 
11799 


377422 
210 
11799 






20 


40 
40 
1 


2663 
78 
98 


35658 
114 
1392 


173394 
148 
6682 


409730 

180 
15095 


587526 
208 
19154 


642827 
234 
19996 


648260 

246 
20438 


648308 

246 
20438 
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